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ABSTRACT 

NGC 6823 is a young open cluster that lies at a distance of ^2 kpc in the Vulpec- 
ula OBI association. Previous studies using CCD photometry and spectroscopy have 
identified a Trapezium system of bright O- and B-type stars at its center, along with 
several massive O- B- and A-type stars in the cluster. We present optical V7?/and near- 
infrared JHK photometric observations, complemented with 6piteer/IRAC archival 
data, with an aim to identify the young low-mass population and the disk candidates 
in this region. Our survey reaches down to /■~22 mag and Kg^l^ mag. There is sig- 
nificant differential reddening within the cluster. We find a bimodal distribution for 
Ay, with a peak at ~3 mag and a broader peak at ~fO mag. We have classified the 
sources based on the [4.5]- [8] color, which is least affected by extinction. We find a 
^^20% fraction of Class I/Class fl young stellar objects (YSOs) in the cluster, while a 
large 80% fraction of the sources have a Class III classification. We have made use of 
the INT Photometric Ha Survey (IPHAS) in order to probe the strength in Ha emis- 
sion for this large population of Class III sources. Nearly all of the Class III objects 
have photospheric (r'-Ha) colors, implying an absence of Ha in emission. This large 
population of Class HI sources is thus likely the extincted field star population rather 
than the diskless YSOs in the cluster. There is a higher concentration of the Class 
I/II systems in the eastern region of the cluster and close to the central Trapezium. 
The western part of the cluster mostly contains Class Ill/field stars and seems devoid 
of disk sources. We find evidence of a pre-main sequence population in NGC 6823, 
in addition to an upper main-sequence population. The pre-main sequence popula- 
tion mainly consists of young disk sources with ages between ^-^1-5 Myr, and at lower 
masses of ~0. 1-0.4 Mq. There may be a possible mass dependent age spread in the 
cluster, with the older stars being more massive than the younger ones. The presence 
of young disk sources in NGC 6823 indicates similar star formation properties in the 
outer regions of the Galaxy as observed for young clusters in the solar neighborhood. 

Key words: Stars: pre-main sequence - Stars: low-mass - circumstellar matter - 
open clusters and associations: individual: NGC 6823 



1 INTRODUCTION 

The open cluster NGC 6823 is considered to be the core of 
the Vulpecula OBI association. The cluster is surrounded 
by the reflection nebula NGC 6820, and is situated in a 
bright H II region, known as Sharpless 86 (Sharpless 1959), 
with associated dark clouds and pillar-like structures (e.g., 
Turner 1979; Chapin et al. 2008; Billot et al. 2010). The 
cluster has a Trapezium system, BD-|-22°3782, of bright O- 
and B-type stars at its center. Erickson (1971) conducted a 



proper-motion membership study of the cluster in a 30' x 30' 
region centered on the central Trapezium system, and deter- 
mined membership probabilities for 92 stars brighter than 
'|/=13 magnitude. Erickson (1971) noted the presence of an 
inner and outer cluster region, such that over one-half of 
the members were found to lie in a region extending to dis- 
tances larger than 4' from the cluster center. Stone (1988) 
conducted a follow-up [/BF photometric study using pho- 
tographic plates of the regions studied by Erickson (1971). 
The H-R diagrams for both the inner and outer cluster re- 
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gions showed a well-defined upper main-sequence, with an 
estimated age of ~2Myr. In addition, Stone (1988) noted 
a pre-main sequence population with ages between 0.2 and 
0.5 Myr in the outer region of the cluster, extending out 
to 15' from the cluster center. Based on the observed age 
dispersion, Stone suggested that NGC 6823 may have expe- 
rienced a recent burst of stellar formation, possibly caused 
by a supernova explosion of one of the massive O-stars at 
the center, with a second generation of stars being formed 
in the outer region some ~0.5 Myr ago. The presence of a 
pre-main sequence population in NGC 6823 has also been 
reported by Guetter (1992), Massey et al. (1995), Pigulski 
et al. (2000) and more recently by Bica et al. (2008), with 
an age spread estimated to be as large as ~2-ll Myr for the 
cluster stars. There have been a few spectroscopic surveys 
of NGC 6823. The earliest type classified yet is 07, and the 
latest is K4 (Turner 1979; Massey et al. 1995; Shi & Hu 
1999). 

The cluster is located at a distance of 2-2.5 kpc (e.g.. 
Turner 1979; Guetter 1992; Massey et al. 1995). The redden- 
ing across the cluster is found to be large and non-uniform. 
The E(B-V) color excesses range between 0.7 and 1.1 mag, 
with a mean value of ~0.9 (e.g., Erickson 1971; Guetter 
1992; Massey et al. 1995; Pigulski et al. 2000). The extinc- 
tion is highest in the eastern part of the central Trapezium, 
in the direction of the reflection nebula NGC 6820 (e.g., 
Pigulski et al. 2000). The reddening for the eleven central 
Trapezium stars is found to be very uniform, with E(B- V) — 
0.84±0.01 mag (Guetter 1992; Pigulski et al. 2000). There is 
evidence of foreground reddening toward NGC 6823, and a 
high contamination by foreground field stars (e.g., Guetter 
1992; Pigulski et al. 2000; Bica et al. 2008). Guetter (1992) 
estimated that the foreground reddening toward the cluster 
is at least 0.6 mag. The presence of a nearby interstellar 
cloud in the direction of NGC 6823 has been confirmed by 
Fresneau & Monier (1999), and Pigulski et al. (2000) suggest 
that more than half of the total absorption towards the clus- 
ter comes from this nearby absorbing matter, at a distance 
0.2-0.5 kpc. 

Chapin et al. (2008) reported the detection of three 
pillar-like structures, VulP-12-13-14, located north-east of 
NGC 6823. Such pillars have been associated with recent 
episodes of star formation (e.g., Chen et al. 2006). The 
close proximity of the supernova remnant SNR G59.5-)-0.1 
to NGC 6823 suggests a possible case of triggered star for- 
mation in this cluster (e.g.. Billot et al. 2010). This is thus 
a young cluster that has likely seen a recent burst of star 
formation. Most of the past surveys in this region have been 
conducted at the wavelengths, and have targeted the 

high mass population. We present deep optical VRI and 
near-infrared JHK photometric observations of NGC 6823, 
complemented with Spitzer/IRAC observations, with an aim 
to identify the low mass population and the disk candidates 
in this region. While searches for young very low mass stars 
have been conducted in several nearby star-forming regions, 
more such deep surveys need to be conducted in distant mas- 
sive clusters. Our observations are presented in Section ij2] 
Section §3.1| discusses the extent of reddening in this region. 
The disk population in NGC 6823 is discussed in Sections 
§3.2| and §3.3| Section §3. 4| presents the spatial distribution 
of the young stellar objects (YSOs) detected, and in Sec- 



tion ij4] we investigate for any observational diagnostics of 
triggered star formation in NGC 6823. 



2 OBSERVATIONS AND DATA REDUCTION 

The optical V, R, and / band images were obtained using 
the Prime Focus camera (WHT/WFC detector) mounted 
on 4m William Herschel Telescope (WHT) in La Palma, 
Canary Islands, Spain. The camera contains two EEV-42- 
80 thinned and AR coated CCDs butted along their long 
axis to provide a 4K x 4K pixel mosaic. The pixel scale for 
this detector is 0.24 arcsec/pixel with a field of view of 16.4 
X 16.4 arcmin. Observations were performed in May 2005, 
obtaining one 3s exposure and one 600s exposure in each of 
the bands. The 3s exposure was to obtain photometry of the 
bright objects which would saturate in the 600s exposures. 
The near-infrared (NIR) J, H, band images were ob- 
tained using the Infrared Side Port Imager (ISPI) mounted 
on CTIO 4m Blanco Telescope in Cerro Tololo, Chile. ISPI 
is a 2K X 2K HgCdTe HAWAII-2 array with a pixel scale 0.3 
arcsec/pixel in the near infrared with a 10.25 x 10.25 arcmin 
field of view. Observations were performed in March 2007. 

The optical and the NIR data were reduced using Image 
Reduction and Analysis Facility (IRAF) system routines. 
For the NIR images, a single stacked image in each band (J, 
H, Ks) was obtained by combining several short-exposure 
images. We conducted Point Spread Function (PSF) pho- 
tometry on the bias subtracted and flat fileded images., us- 
ing various routines under the IRAF/ DAOPHOT package 
(Stetson 1987). The routine daofind is first applied to de- 
tect all the stellar sources. The computation of the PSF to 
be fitted was performed by the routine psf. The PSF was 
modeled to be the sum of an analytic bivariate Gaussian 
function and empirical corrections from the best Gaussian 
of the true observed brightness values within the average 
profile of several stars in the image. PSF photometry was 
performed with the routine allstar, which after classifying 
the stars in groups, compiles a catalog of the most likely 
candidate stars, based on their PSF fitting and the physical 
conditions, and subtracts them from the original image. The 
NIR absolute photometry was calibrated using the 2MASS 
photometric system. The optical photometry was calibrated 
using sources from CADC photometric standard list for the 
NGC 6823 field. Astrometric solutions were obtained for the 
optical and NIR frames using imwcs command in WCSTools 
(Mink 1999). We then created an optical -I- NIR catalog using 
gcompare, a tool from the ELIXIR system. Separate cata- 
logs were created for the bright and the faint sources, and 
any duplicates were removed. We then cross matched these 
catalogs with data from the Spitzer GLIMPSE survey. Our 
final catalog consists of 627 sources, with good photomet- 
ric detections (cr <0.1 mag) in all four IRAC bands. For 17 
of these sources, we have Spitzer MIPSGAL 24ii.m matches 
obtained from Billot et al. (2010). None of the sources were 
detected at 70iJ.m. Billot et al. estimate an uncertainty of 
<0.2 mag for the 24|i.m photometry. The optical, NIR and 
Spitzer/IRAC photometry for these sources is listed in Ta- 
ble [l] and the MIPS 24nm photometry is listed in Table [2] 

Figure [T] gives a spatial view of NGC 6823. This is an 
IRAC composite image of 3.6(J.m (blue), 4.5vim (green), and 
8vim (red) observations. The survey area from our work is 
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marked by white lines, the central Trapezium is marked in 
black. Taylor et al. (1992) detected a supernova remnant, 
SNR G59. 5+0.1, in the direction of the Vulpecula OBI as- 
sociation. The location of this SNR and the three pillars, 
VulP-12-13-U, detected by Chapin et al. (2008), are also 
marked. 



3 ANALYSIS 
3.1 Extinction 

Figure [2] shows the NIR color-color diagram (ccd) for the 
NGC 6823 cluster. The thick black and red curves indicate 
the zero-age main sequence for dwarfs and giants, respec- 
tively (covering spectral types 0-M9). The zero-age main 
sequence was assembled from the work of Bessell & Brett 
(1988), Leggett (1992), Leggett, AUard & Hauschildt (1998), 
and Kirpatrick et al. (2000). The blue line marks the clas- 
sical T Tauri star (CTTS) locus from Meyer et al. (1997). 
The dashed lines are parallel to the reddening vector, and 
intersect the dwarf main-sequence at spectral types of MO 
and M9. The NIR ccd shows a highly extincted population of 
early-type main-sequence dwarfs concentrated near (H—Ks) 
~ 0.7, {J — H) ~ 1.5. These stars form a sequence parallel to 
the reddening vector, and comparing with the length of the 
Av vector, we can infer that these are located at Av 5S 10. 
There is another group of stars that lies close to the CTTS 
locus near (H - K,) ~ 0.6, {J ~ H) 1. This is hkely the 
young disk population in the cluster. We also find a more 
dispersed population that have modest Ks-hand excesses, 
with (H — Ks) J? 1. This 'red' population lies to the right 
of the principal sequence of reddened field stars, and has a 
different slope than the reddening vector shown. 

The median (H — Ks) color for YSOs is only about 
~0.2 magnitude (e.g., Wilking et al. 1999), which makes 
it difficult to distinguish young stars with modest Ka-hand 
excesses from the extincted field stars. Such a distinction 
can be made more clearly if a combination of NIR and mid- 
infrared (MIR) colors is used. In order to estimate the extent 
of reddening in NGC 6823, we have considered the Rayleigh- 
Jeans Color Excess (RJCE) method described in Majewski 
et al. (2011). This method uses a combination of NIR and 
MIR measured colors, and extinction to a star is determined 
by dereddening the observed colors along color excess vec- 
tors to the intrinsic stellar locus. Majewski et al. have based 
the RJCE method on the {H - [4.5]) color, since the spread 
in the intrinsic {H - [4.5] )o color is found to be very small (sC 
0.4 mag) over a wide range of B through M spectral types. 
Due to this near constancy, the interstellar reddening can 
be estimated using only the observed {H - [4.5]) color. The 
RJCE conversion relation is, Ak = 0.918 {H - [4.5] - 0.08), 
where Majewski et al. (2011) have considered an intrinsic 
{H - [4.5])o color of 0.08, and the reddening law from Inde- 
betouw et al. (2007) has been used to derive the conversion 
relation. 

Figure ^ shows the Av histogram for NGC 6823 {Ak 
= 0.112 Av; Rieke & Lebofsky 1985). We find a bimodal 
shape for the Av distribution, with a peak at ylv~3 mag 
and a broader peak at Av^W mag. To demonstrate the bi- 
modality, we applied a simple goodness-of-fit test. Our 
null hypothesis is that a unimodal distribution cannot ad- 
equately represent the observed Av distribution. We fitted 



the Av histogram by a Gaussian distribution, using the MP- 
FITPEAK program in IDL. The Gaussian fit is shown in 
Fig. [3^ (bottom panel). The peak value of the Gaussian is 
42.64 and the Gaussian sigma is 7.07. The number of de- 
grees of freedom (dof) is computed as the difference of the 
number of elements and the number of free parameters used 
for the fit. The Gaussian peak value (height), sigma (width), 
and the centroid were used as the free parameters for fitting. 
The resulting (summed squared residuals) value of the fit 
is found to be 1745.06, and the reduced-x^ value (x^/dof) 
is 32.3. For an ideal fit, we would expect the reduced-x'^ 
value to be less than or near one. The large value obtained 
indicates that the fit has missed too many data points and 
is a poor fit. The associated x^ probability, Q, which esti- 
mates the significance of the x^ test, is found to be 0.0. A 
high value of Q indicates that the resulting x^ value will be 
different if the test is repeated again. If the observed Av dis- 
tribution were unimodal with a single broad peak observed 
around Av~10 mag, we would have obtained a good fit with 
the Gaussian distribution and the resulting reduced-x^ value 
would be close to one. Thus our null hypothesis is correct 
and the bimodality in the Av distribution is likely real. 

A large fraction of the sources are part of the more ex- 
tincted population with Av>5 mag. Fig. [sJj compares the 
spatial distribution of these two groups of sources, with Av 
^ 5 mag (blue points) and Av >5 mag (green points). The 
central Trapezium in NGC 6823 consists of 11 bright O- 
and B-type stars, which show nearly uniform reddening of 
Av = 2.57±0.031 mag (Guetter 1992; Pigulski et al. 2000). 
The Trapezium stars are thus part of the smaller group of 
sources with a peak at Av~3 mag. Overall, we do not find 
any prominent high- or low- extinction regions in the clus- 
ter. Pigulski et al. (2000) had noted higher extinction in the 
eastern part of the central Trapezium, and had explained it 
by the presence of the reflection nebula NGC 6820 in that 
direction. In Fig. [sj; (top panel), we have compared the ex- 
tent of extinction in the eastern and western parts of the 
cluster. We have plotted here the separation in RA for each 
source from the central Trapezium. We have considered the 
star E81 Sh-d [identification numbers from Erickson (1971) 
and Sharpless (1959)] located at (q, S) = (295.8°, 23.3°) to 
be the center of the Trapezium system. This is a B1.5V star, 
and has a 98% cluster membership probability from Erick- 
son (1971). In Fig. [3]:, positive values for ARA imply sources 
located to the east of the central Trapezium. A wider range 
in Av between and 50 mag is observed in the eastern 
region, while sources in the western region show two con- 
fined distributions at ~0 and ~10 mag. The bottom panel 
in Fig. |3]; compares the cumulative distribution functions, 
obtained from a Kolmogorov-Smirnov (K-S) statistical test, 
for the high-extincted {Av>5 mag) and the low-extincted 
(j4v^5 mag) sources located to the east (solid line) and 
to the west (dashed line) of the central Trapezium. The 'X' 
variable in these distributions implies the values for Av ■ For 
the low-extincted sources {bottom panel, left), the east and 
the west distributions have similar shapes. According to the 
K-S test, the maximum difference between these cumula- 
tive distributions is 0.11, and the corresponding probability 
that the two distributions are similar is P=0.7. In compar- 
ison to this, the cumulative fraction of the high-extincted 
sources with Av between ~10 and 20 mag is slightly higher 
in the western part of the cluster than the east. This can be 
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Figure 1. A spatial view of NGC 6823 and the Vulpecula association. This is an IRAC composite image of 3.6iJ.m (blue), 4.5|im (green), 
and 8|^m (red) observations. The color intensity bar in units of MJy sr~^ is shown at the bottom of the image. North is up, east is to the 
left. The region in NGC 6823 that we have surveyed is marked by a white rectangle, and is centered on the Trapezium (marked in black). 
In the north of the Trapezium is the supernova remnant SNR G59. 5+0.1. Also marked are the three pillar-like structures, VulP-12-13-14 
(green diamonds). 



seen from the bottom, right panel in Fig.jSf;. The difference 
between the two cumulative distributions is not significant; 
according to the K-S test, the maximum difference is 0.24, 
and the probability P that the two distributions are simi- 
lar is 0.6. There may be a slightly higher concentration of 
the high-extincted sources located west of the Trapezium, 
but no significant spatial concentration of the high- or low- 
extincted sources is evident. Nevertheless, there is signifi- 
cant differential reddening within the cluster, and we can 
infer that the cluster lies behind at least ~3 magnitudes of 
visual extinction. 

The main uncertainties associated with the RJCE 
method are due to the spread in the intrinsic stellar colors, 
the uncertainties in the extinction law used, and the uncer- 
tainties in the photometric measurements. Majewski et al. 
(2011) estimate an uncertainty on Ak of ^0.1 mag. We note 
that the intrinsic {H - [4.5])o color of 0.08 used in their re- 
lation is for F, G and K stars, and the intrinsic color could 
be as large as 0.4 mag if the later type stars are taken into 
account. This can lead to a difference of ~3 magnitude in 
the Av values, which is small in the context of Fig. [2^. The 
intrinsic variability of a star makes it difficult to accurately 
estimate the extinction towards each star individually, but 
such methods are reliable in estimating the extent of red- 
dening over the full survey area studied. A commonly used 
method of estimating Av is by dereddening each source in 



the {H-Ks)/{J-H) ccd to the CTTS locus (e.g., Luhman 
et al. 2006). This method, however, could lead to large un- 
certainties if the sample consists of field stars and diskless 
weak-line T Tauri stars (WTTS). Also, the spread in the 
intrinsic NIR colors that include the J-band photometry is 
larger than the NIR-I-MIR color combinations (e.g., Majew- 
ski et al. 2011). The RJCE provides a more reliable method 
for estimating extinction, compared to previous dereddening 
methods based on NIR colors only. 



3.2 IRAC Source Classification 

We have used the empirically-derived methods described in 
Gutermuth et al. (2008) to classify the sources that show 
excess emission in the mid-infrared IRAC bands. There are 
a few well-known classification schemes that use the 2-8jim 
slope to classify YSOs into the different evolutionary classes 
(e.g., Allen et al. 2004; Lada et al. 2006). These schemes 
include the 3.6[i.m photometry, which is strongly affected by 
reddening from dust. From the extinction law by Indebe- 
touw et al. (2005) and Flaherty et al. (2007), the extinction 
is nearly constant at 0.43A_ff in the 4.5jim through 8jim 
bands, but rises to 0.56A_fc in the 3.6|J.m band. Gutermuth 
et al. (2005) have made use of the flattening of the extinc- 
tion law in the IRAC 4.5 to 8iJ.m bands, and have developed 
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Figure 2. The near-infrared ccd for NGC 6823. The thick black and red curves indicate the zero-age main sequence for dwarfs and giants, 
respectively. The thin black line marks the CTTS locus from Meyer et al. (1997). The reddening vector shown is from the extinction 
law of Reike & Lebofsky (1985). The black dashed lines are parallel to the reddening vector, and intersect the dwarf main-sequence at 
spectral types of MO (left line) and M9 (right line). Typical l-cr errors in (J-H) and (H-Ks) are 0.06 mag and 0.02 mag, respectively. 



a classification scheme based on the [4.5]-[5.8] color to dis- 
tinguish between Class I and Class II YSOs. This color is 
much less affected by dust extinction than any [3.6]-based 
color. 

The first phase in their classification method is to re- 
move the various extragalactic contaminants. After applying 
the respective color cuts given in Gutermuth et al. (2008) 
for identifying possible extragalactic contaminants, we have 
found 137 sources which have colors consistent with being 
PAH, AGN, or shock emission sources. After removing these 
contaminants, the remaining 490 sources were classified into 
Class I or Class II by applying the appropriate color con- 
straints defined in Gutermuth et al. (2008). In total, we find 
94 Class II and 2 Class I systems. The remaining sources 
have colors consistent with photospheric emission, and were 
classified as Class III sources. We note that objects lying 
in this Class III category can be field stars and/or diskless 
YSOs. They can also be transition disks with significant in- 
ner disk clearing or large inner holes in the disk. Such sources 
exhibit photospheric emission in the IRAC bands, but could 
flare up at longer wavelengths. In the discussion that con- 
tinues, we will refer to the Class III sources as Class Ill/field 
stars. In total, these Class Ill/field stars are 394 in number. 
We thus have a ~20% fraction of Class I and Class II YSOs, 
while about 80% are Class Ill/field stars. 

Figure |4|i shows the {H — Ks) versus [4.5]-[8] ccd for 
the 490 sources. The YSOs and the diskless sources show a 
similar wide range in the NIR colors, but the IRAC color 
provides a good separation between the photospheric and 
the disk population at [4.5]-[8]~0.4. There are a few Class 
Ill/field stars with IRAC colors redder than this value, and 
are more mixed with the Class II sources. These could be 
the Class III YSOs in the cluster, that show slightly redder 
IRAC colors than the main-sequence population. The IRAC 



ccd in Fig. |4j3 shows a nice separation between the different 
evolutionary classes. The protostars are the reddest in the 
[3.6]-[5.8] color and lie in the top right corner. The Class 
II systems can be identified at [4.5]-[8]^0.5, [3.6]-[5.8]^0.4, 
while the Class Ill/field stars form a locus near [4.5]-[8] and 
[3.6]-[5.8]~0.2-0.3. 

In Fig. |4|;, the inner box indicates the IRAC colors ob- 
served for the T Tauri stars in Taurus, from the work of 
Hartmann et al. (2005). The median IRAC colors in Taurus 
are [3.6]-[4.5] = 0.4, [4.5]-[5.8] = 0.5, [5.8]-[8] = 0.8. For the 
Class II sources in NGC 6823, the median colors are [3.6]- 
[4.5] = 0.5, [4.5]-[5.8] = 0.3, [5.8]-[8] = 0.7. Most NGC 6823 
Class II YSOs show stronger excess in the [3.6]-[4.5] color 
compared to Taurus, which can be explained by the higher 
extinction observed in this cluster. Taurus, in comparison, 
is relatively free of extinction {Av^2 mag; e.g., Luhman et 
al. 2006). 

We have MIPSGAL matches from Billot et al. (2010) 
at 24fi.m for 15 Class II YSOs, and one Class Ill/field star. 
Fig |4jl shows the IRAC-MIPS ccd for these 16 disks. The 
Class II YSOs show larger excesses in the [8]- [24] color 
than [3.6]-[5.8], which indicates some flaring in the outer 
disks for these systems. The [8]-[24] color for the one Class 
Ill/field star is consistent with photospheric emission. Fig. [5] 
shows the spectral energy distributions (SEDs) for these 
16 sources. The top left SED is for the one Class Ill/field 
star (194312.47-1-231817.03). For most Class II sources, the 
shape of the SED has either flattened between 8 and 24/im, 
or shows some flaring between these two points. The YSO 
194308.33+231913.01 shows strong flaring at 24nm. The 
MIPSGAL survey has a 5-(j rms point source sensitivity of 
1.3 mjy at 24/im. Thus only the brightest disks could be 
detected by this survey. Deeper observations should result 
in a larger number of detections. 
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Figure 3. (a): (Top) The Ay histogram for NGC 6823 sources, obtained using the RJCE method. A bimodal distribution can be seen, 
with peaks at ~3 and ~10 mag. Bottom panel shows a Gaussian fit to the observed Ay distribution. More details are provided in the 
text. 



3.3 Ha Emission Strength 

One way to distinguish between Class III YSOs and field 
stars is the strength in the Ha emission. We matched our 
NGC 6823 catalog with the INT Photometric Ha Survey 
(IPHAS) catalog (Drew et al. 2005). We were able to obtain 
IPHAS matches for 213 sources. Fig. |6] plots the IPHAS 
(r' — i') versus (r'-Ha) colors. Also overplotted is the grid 
of simulated IPHAS colors from Barentsen et al. (2011) for 
spectral types of F0-M6. These simulated colors track the 
increase in the Ha emission for a given spectral type. We 
have not applied any extinction correction to the observed 
colors or the simulated grid. The (r' — i') color thus rep- 
resents a combination of the interstellar reddening and the 
intrinsic color of the source. The (r'-Ha) color gets redder 
with increasing emission or larger equivalent width (EW) in 
the Ha spectral line. The reddening vector shown in Fig. [6] 



is from the optical-IR extinction law of Howarth (1983) and 
Schlegel et al. (1998). 

Nearly all Class Ill/field stars lie below the OA EW 
track, implying an absence of Ha in emission and photo- 
spheric (r'-Ha) colors for these sources. We find a concen- 
tration of these photospheric sources at (r' — i') of ~0.3-l, 
and these are likely the FGK type field stars. There is also 
a dispersed population at (r' — i') ^2 that lies below the OA 
track. The red (r' — i') colors for these sources are likely af- 
fected by high extinction, as these lie in cluster regions where 



Av>7 mag (Section f3.4l. We expect the diskless sources 



that lie below the OA track to be main-sequence stars. 

Young disk sources are expected to be CTTS with on- 
going accretion in the disk. A boundary of Ha EW=10 A 
is traditionally used to make a distinction between CTTS 
and the non-accreting WTTS (e.g., Barentsen et al. 2011). 



Young Stellar Objects m NGC 6823 7 




-OAl aS2 1.5 2,4 3_3 4.3 5 3 6.1 7.1 

Figure 3. (b): A comparison of the spatial distribution of the two groups of sources with Ay ?J 5 mag (blue points) and Ay >5 mag 
(green points). The central Trapezium is marked by black lines. 



Most of the Class II systems in Figjojlie above the lOA track, 
and can be classified as accreting CTTS systems. We also 
have about 20 Class II sources that lie below this bound- 
ary, which indicates the presence of passive or non-accreting 
disks around these sources. Emission in Ha arises from the 
hot spots on the stellar surface, where the accreting ma- 
terial from the circumstellar disk strikes onto the star at 
free-fall velocities (e.g., Hartmann et al. 2005). The origin 
of Ha emission in disk bearing sources is from the inner disk 
regions close to the central source (<1 AU for a typical T 
Tauri disk) . The Ha emission line is thus a gas diagnostic of 
the warm inner part of the disk and is not indicative of the 
total gas content in the disk. The few Class II cases with an 
absence of Ha emission may be more evolved systems with 
inner disk clearings. If the disk material has dissipated in 
the inner regions then we would not expect emission in the 
Ha line. Fig. [6] indicates that selecting YSOs based on the 
strength in Ha emission could result in a rejection of young 
disk sources with inner disk clearings. 

Among the Class III sources, there are 8 that have Ha 
EW larger than lOA and may be classified as CTTS sys- 
tems. Emission in Ha, however, could also arise from an 
active chromosphere. Chromospheric activity in field stars 
is found to show a rise towards later types, with a peak 
at spectral type of M6-M7 (e.g., Riaz et al. 2006a). The 
8 Class HI sources with strong Ha emission have (r' — i') 
colors similar to the later spectral types of M4-M6. These 
could be chromospherically active field stars, rather than 
young CTTS objects. 

Wc could select all CTTS sources as candidate clus- 



ter members, and reject the sources without Ha emission, 
particularly the Class Ill/field stars with photospheric (r'- 
Ha) colors as background/foreground field stars. However, 
some of these photospheric sources may be diskless YSOs. 
We plan to conduct an X-ray imaging survey of NGC 6823 
which will be important in identifying the genuine YSOs 
among the Class Ill/field star population. The IPHAS pho- 
tometry for the 213 sources is listed in Table |3j with the 
possible accretors marked by an asterisk. 

3.4 Spatial Distribution 

Figure [T] shows the spatial distribution for the different SED 
classes in NGC 6823. We find a higher concentration of Class 
II systems in the eastern region of the cluster and close to 
the central Trapezium. The two protostellar sources are lo- 
cated in the northern part, with one protostar located close 
to the VulP12 pillar in the north-east. The western part of 
the cluster mostly contains Class Ill/field stars. The differ- 
ences in the spatial distribution of the Class I/II and diskless 
population can be noted more clearly in Fig.jS] Most Class II 
systems lie within 100" of the central Trapezium (top panel, 
left). Another peak is seen at a distance of ~200" east of the 
Trapezium stars. These are two groups of Class II sources, 
one located near the VulP12pA\&r observed in the north-east 
and another more dispersed group located in the south-east. 
The Class Ill/field stars are more evenly distributed in the 
cluster, though a higher number of these sources are found 
in the western region (Fig.jsj top panel, right). Fig. [8]( bottom 
panel) shows the Av distribution for the Class I/II and the 
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Figure 3. (c): (Top) The separation in RA from the central Trapezium vs. Ay- ARA = implies the location of the central Trapezium, 
positive values for ARA imply location east of the Trapezium. Bottom panel compares the cumulative distribution functions, obtained 
from a K-S test, for the low-extincted (left) and the high-extincted (right) sources located to the east (solid line) and to the west (dashed 
line) of the central Trapezium. Here, 'X' implies the values for Ay. More details are provided in the text. 



Class Ill/field stars. The disk population is mostly located in 
the high extinction regions, where Av is ~10-20 magnitude. 
The diskless sources show a bimodal Av distribution, similar 
to that observed for the whole cluster (Fig. [2|, which indi- 
cates their large spread over the whole region. The western 
region of NGC 6823 thus seems more devoid of disk sources. 
This western part could be more evolved, so that the disks 
have all been dissipated, or could have a high contamina- 
tion of foreground/background field stars. Fig. [o] shows the 
spatial distribution of the possible accretors in NGC 6823 
(Hof EW > 10 A). Most of the accretors among the Class III 
sources are located in the southern region, while the Class 
II CTTS sources are more dispersed, but are located mainly 
near the Trapezium or the eastern part of the cluster. The 
main conclusion that can be made from Figs. [7] and [9] is that 
star formation is not concentrated near the center in NGC 
6823, and has spread to the outer regions, mainly towards 
the east of the central Trapezium. 



3.5 Age and Mass Estimates 

Figure[lO]shows the V versus (V-I) color-magnitude diagram 
(cmd) for NGC 6823. We have considered the pre-main se- 



10 



quence models from Siess et al^(2000; SOO) (Fig.[lOji), Palla 
& Stabler (1999; PS99)_(Fig 
(1998; DAM98) (Fig 



IOJd) and D'Antona & Mazitelh 
to estimate the ages and masses 
for our sample. The evolutionary tracks and isochrones are 
for solar metallicity, and have been scaled to a distance of 
2 kpc. A first analysis of the cmds from all three models 
indicates the possible presence of two different populations 
in the cluster. One is a sequence of older stars located along 
the zero-age main sequence. These sources are also the more 
massive in the cluster, with M* > 0.4 Mq. A majority of 
these objects are Class Ill/field stars, with photospheric (r'- 
Ha) colors, and {r'-i') colors consistent with FGK type stars 
(Fig.jef. The second group consists of younger sources with 
ages of <10 Myr, and at lower masses of ~0. 1-0.4 Mq. A 
large fraction (~75%) of these are Class II systems. There is 
also a smaller group of mainly Class Ill/field stars that lie 
at the faint end, with V J? 20 and (V-I) ^ 4.5 mag. A few 
of these could be the candidate young brown dwarfs in the 
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Figure 4. IRAC ccds for NGC 6823 sources. Red points denote Class I systems, blue are Class II and grey are Class Ill/field stars. 
The reddening vectors are plotted using Indebetouw et al. (2005) extinction law. The boxed area in the bottom panel (left) denotes the 
median Taurus colors. 



cluster. In Fig. |10[ we have marked with circles the sources 
that may be possible accretors. Most of the Ha emitters are 
among the younger group of sources, while the Ha emitters 
among the older group are all Class H systems. The zero- 
age main sequence was assembled from the work of Bessell & 
Brett (1988), Leggett (1992), Leggett, AUard & Hauschildt 
(1998), and Kirpatrick et al. (2000). 

Figure [Tl] shows the resulting age and mass distribu- 
tions using the three different models. We have labelled the 
oldest age bin as ">10 Myr". Though we expect most of 
these sources to be the main sequence population in the 
cluster, this older bin also includes about 25 Class H sys- 
tems that lie between the 10 Myr isochrone and the zero- 
age main sequence. Being disk sources, these are expected 
to be younger than the main-sequence population. These 
older disk sources are likely the secondary debris disk sys- 
tems, with a regeneration of warm dust in the inner disk 
regions (e.g., Riaz et al. 2006b). The age distributions using 
the PS99 and DAM98 models indicate a strong peak for the 
older sources, along with a tail of objects at younger ages of 
<10 Myr. The shape of the age distribution using the SOO 
model is different, with two comparable peaks at 1-5 Myr 
and >10 Myr. There is a wider range in the V magnitudes 
between the 1 and 5 Myr isochrones for the SOO model com- 
pared to the other two cases (Fig. 10 1), which would explain 



the two comparable peaks observed for the SOO age distribu- 
tion. From all three models, most of the young sources seem 



to be at ages of ~0.1-5 Myr, and the least number of objects 
is found in the 5-10 Myr age bin. For the PS99 model, the 
0.1 Myr isochrone is ~2 magnitudes brighter than the 1 Myr 
isochrone at a given (V-I) color, while the diflerence is less 
than one magnitude for the SOO or the DAM98 models. This 
is particularly for objects redder than (V-I)~3. This results 
in a greater number of objects in the 0.1-lMyr bin for the 
PS99 model compared to the other two models. We note 
that for all three models, the <0.1 Myr age bin may not be 
reliable. This bin mainly contains the group of faint sources 
at V ~ 22 and (V-I) 4.5 mag, most of which are Class 
Ill/field stars. It is unusual for a source which is presumably 
at an age of <0.1 Myr to be devoid of circumstellar material. 
These red objects could be extincted field stars rather than 
very young low-mass stars. Also, for the DAM98 model, the 
isochrones are not available for (V-I) bluer than ~1.5. We 
have roughly extrapolated the 1, 5 and 10 Myr isochrones 
from this model to estimate the ages for these bluer sources, 
which are likely to be older, massive objects. Thus the 1-5 
Myr and >10 Myr age bins for the DAM98 model could have 
fewer number of sources than plotted in Fig. |11[ 

Other than these three models, we have also considered 
the evolutionary models by Baraffe et al. (1998; BC98). The 
cmd and the resulting age distribution from this model are 
shown in Fig. |12[ The BC98 model does not provide a 0.1 
Myr isochrone. As an be seen, nearly all of the young sources 
lie above the 1 Myr isochrone. The resulting age distribution 
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Figure 5. SEDs for the 16 sources with 24/im MIPSGAL matches. 



shows two strong peaks, one for < 1 Myr age bin, and the 
other for the older hkely main-sequence population. It is very 
difficult to estimate masses using this model, particularly 
for the younger sources extending from (V-I) 3. We have 
therefore not attempted to build the mass distribution using 
the BC98 model. 

Despite the uncertainties and the differences among the 
models, we find a pre-main sequence population in NGC 
6823, in addition to an upper main-sequence population. 
The age distributions for NGC 6823 indicate recent star for- 
mation activity in the last ~l-5 Myr, and there may be 
candidate cluster members at ages as young as 0.1 Myr. The 
star formation history of NGC 6823 seems similar to that ob- 
served in young star-forming regions such as the Orion Neb- 
ula Cluster, Taurus-Auriga, Lupus and p Ophiuchi, which 
also show an older population of ~10^ yr, along with more 
recent star formation activity in the last 1-2 Myr (e.g., Palla 
& Gaali 1997; Martin et al. 1998; Muench et al. 2002). The 
similarities with NGC 6823 suggest that star formation may 
proceed in more than one epoch within a wide range of en- 
vironments. 

The stellar mass tracks from these different evolution- 
ary models show a more prominent difference than the 
isochrones (Fig. 



10 1. Based on the SOO model, all of the 



sources in the cluster have masses higher than 0.1 A/q 
(Fig. [Tok), while using the DAM98 model results in a few 



parison, the lowest mass track from the PS99 model is for 
O.IA/q, and it lies at a bluer color than the 0.1 Mq track 



from DAM98 or SOO model (Fig. 10 a). Thus the PS99 model 



results in a larger number of sources with masses less than 
O.lAfQ. Also notable are the nearly constant (V-I) colors at 
a given mass for the PS99 case, while the mass tracks for 
M,^0.2 Mq for the DAM98 model get redder with increas- 
ing age (Fig. 10d,c). The resulting mass distributions from 



these models thus have very different shapes (Fig. 11 right 
panel). The mass distribution from the SOO model shows 
a prominent peak for massive sources at M*>1Mq, with 
a nearly flat distribution for lower masses and then a rise 
again for the 0.1-0.3 Mq mass bins. The PS99 model (mid- 
dle, right panel) also shows peaks at the two ends of the 
mass distribution, with a flat distribution for masses of 0.2- 
1 Mq. The distribution based on the DAM98 model also 
shows a peak for the very low-mass and the very high-mass 
sources, but there is almost a break in the distribution at 
O.3-O.4M0. The range in the (V-I) colors between the 0.3 
and 0.4 Mq tracks from this model is very narrow (Fig. 



10 



which results in a negligible number of sources in this mass 
bin. This is also the case for the 0.8-0.9 AIq bin. On the 
other hand, the range in the (V-I) colors between 0.1 and 
0.2 Mq is much larger for DAM98 than the other two mod- 
els (Fig. 10:), which results in a higher number of sources 



sources having lower masses of ~0.08AfQ (Fig. lOi). In com- 



in this mass bin compared to the other two models. 

We note that there are several inconsistencies in the 
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Figure 6. The (r' — «') versus (r'-Ho) ccd for NGC 6823. Also overplotted is the grid of simulated IPHAS colors from Barentsen et 
al. (2011) for spectral types of F0-M6 and Ho EW of 0-200A. Red points denote Class I systems, blue are Class II and grey are Class 
Ill/field stars. 



mass tracks available. The boundaries between the evolu- 
tionary tracks for masses higher than ~0.4Mq are less clear 
for the older sources at ages >10 Myr (Fig. 10 1. We have 



roughly extrapolated the mass tracks to estimate masses for 
these older sources. We would thus consider the mass bins 
of >0.4 Mq to be unreliable. For the DAM98 model, the 
highest mass track is for 0.9Mq, and so the >0.9A/q bin 
for this model includes all of the sources that lie above this 
track (Fig. 
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On the other hand, the >1 Mq mass bin 
for the PS99 and SOO models includes sources with masses 
of 1-3 Mq and 1-5 Mq, respectively. This makes it diffi- 
cult to directly compare the high-mass peaks observed for 
the three models. It is also difficult to estimate the masses 
for the tail of faint sources observed at V~22 and (V-I) 5s 
4.5. We have included all of these faint objects in the <0.1 
Mq bin, which makes this bin unreliable as it may contain 
fewer sources than plotted. Another inconsistency is that the 
mass tracks from every model are not available at intervals 
of O.IMq, due to which the bin sizes in Fig. 
same. 
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are not the 



Overall, the mass distributions indicate a large spread 
in the stellar masses, with two main groups at ~0.1-0.4 Mq, 
and at >1 Mq. There may be a possible mass-age relation- 
ship in the cluster, or a mass dependent age spread, with 
the older stars being more massive than the younger ones. 
None of these pre-main sequence models, however, can ex- 
plain the small group of faint sources located at (V-I)5s4.5. 
The position of a young star in a cmd is dependent on its 
mass, age, extinction, binarity, its intrinsic variability, and 



the presence of a circumstellar disk. These properties need 
to be taken into account in order to construct the IMF and 
study the star formation history of the cluster. Follow-up 
spectroscopic observations will allow us to confirm cluster 
membership for these potential low-mass YSOs, and to build 
up a more reliable age and mass distribution. 



4 DISCUSSION: A POSSIBLE CASE OF 
TRIGGERED STAR FORMATION 

Taylor et al. (1992) detected a supernova remnant, SNR 
G59. 5-1-0.1, in the direction of the Vulpecula OBI associa- 
tion. This SNR is located at a distance of ~11 pc from the 
central Trapezium in NGC 6823 (Fig. [T}. Taylor et al. de- 
scribed it as a shell-type SNR with a diameter of 15'. Billot 
et al. (2010) noted a high density of YSOs lined up along two 
arcs at the northern and southern rims of SNR G59. 5-1-0.1, 
and suggested that the interaction of the expanding shell of 
SNR G59. 5-1-0.1 and the neutral gas of Sh2-86 could have 
triggered the formation of young stars in the cluster. As 
explained in Preibisch & Zinnecker (2006), close to a su- 
pernova, the shock wave will destroy the surrounding cloud, 
but at large distances of ~20-100 pc, the shock velocities are 
below ~50 km s"'^ and the shock properties are suitable for 
triggering cloud collapse. The proximity of SNR G59. 5+0.1 
to the central Trapezium suggests that the shock waves from 
this explosion may have triggered the star formation activity 
that we witness in the present epoch. An important aspect 
of supernova-triggered star formation is that the age spread 
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Figure 8. Top panel: The separation in RA from tlie central Trapezium for ttie Class I/II (left) and the Class Ill/field stars (right). 
Bottom panel: The Ay distribution for Class I/II (left) and Class Ill/field stars (right). 
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Figure 9. The spatial distribution for accrctors (Ha EW>10Jl). Class II systems are denoted by blue, Class Ill/field stars are denoted 
by green symbols. 



Seiss 




2 4 6 

V-l 



Figure 10. (a): The V vs (V-I) cmds for NGC 6823. The evolutionary tracks and isochrones are from Seiss et al. (2000) models. The 
reddening vector is plotted Reike & Lebofsky (1985) reddening law. Grey symbols indicate Class Ill/field stars, blue symbols denote 
Class II systems. The objects marked with black circles show emission in Ha (EW >10 A). Black thick line is the zero age main-sequence. 



14 Riaz et al. 



PS99 




2 4 6 8 

V-l 

Figure 10. (b): The V vs (V-I) cmds for NGC 6823. The evolutionary models are from Palla & Stahler (1999). Symbols are the same 
as in Fig. |10^ . 
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Figure 10. (c): The V vs (V-I) cmds for NGC 6823. The evolutionary models are from D'Antona & Mazitelh (1998). Symbols are the 
same as in Fig. |10^ . 
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Figure 11. The age {left) and mass [right) distributions for NGC 6823, obtained using the models by Seiss et al. (2000) {top 
Palla & Stahler (1999) {midle panel), and D'Antona & MaziteUi (1998) {bottom panel). The error bars are the l-cr Poisson uncertainties. 



in the new subgroup formed should be small, since the shock 
waves disperse through the cloud on a short timescale. The 
estimated spread in ages for the new generation of stars in 
NGC 6823 is ~l-5 Myr, which is not significantly large. How- 
ever, Billot et al. (2010) estimate an age of ~10'' yr for SNR 
G59. 5+0.1. This event seems too recent to have triggered 
any of the star formation that we see now. 

NGC 6823 is also known to host some 20 spectroscopi- 
cally confirmed OB stars (Fig. [7| , which have cluster mem- 
bership probabilities of >95% (e.g., Shi & Hu 1999). The 
ionizing radiation and winds from massive stars can disperse 
the clouds in the immediate surroundings, and trigger the 
collapse of a molecular cloud which may otherwise not con- 
tract and fragment spontaneously. The classical model for 



triggered star formation from OB stars is the "coUect-and- 
collapse" model (Elmgreen & Lada 1977). In this model, as 
soon as the first massive star is born deep inside a cloud, 
the ionization fronts from the H II regions created by the 
massive star push the cloud from within, forming a cavity. 
The gas and dust in the cloud accumulate to a layer until a 
critical density is reached for gravitational collapse, trigger- 
ing the formation of the next generation of stars. According 
to this model, low-mass stars should be systematically older 
than the massive OB stars, and there should be a large age 
spread, as large as the lifetime of the cloud (e.g., Preibisch & 
Zinnecker 2006). Another observational evidence of trigger- 
ing by this mechanism is the presence of large-scale cavities 
in the association, as observed in the W5 region (W5 West; 




E 80 - 
E 



BC98 



1 

5-10 

Afle (Mvr) 

Figure 12. The V vs (V-I) cmd using Baraffe et al. (1998) models. Bottom panel shows the resulting age distribution. 



Koenig et al. 2008). Such large-scale cavity-like structures 
have not been observed in the Vulpecula OBI association 
(Billot et al. 2010) or in the NGC 6823 region surveyed. 
We do find evidence of the lower mass population to be 
younger than the massive stars in the cluster, but the age 
spread is not as large as the lifetime of the cloud. The coUect- 
and-coUapse model thus does not seem to be the applicable 
mechanism in this cluster. 

Another mechanism of triggered star formation by OB 
stars is called the "radiation-driven implosion" (RDI) pro- 
cess (e.g., Hester & Desch 2005). According to this model, an 
OB star drives an ionization shock front into the surround- 
ing cloud, and cores within the cloud are thus triggered into 
collapse by the shock wave. The RDI model predicts that the 
low-mass stars should be younger than the OB stars, since 



the OB stars initiated their formation. Chen et al. (2006) 
have discussed some more observational diagnostics of this 
model. The remnant cloud will have the shape of an ex- 
tended pillar-like structure, and will extend towards or point 
towards the massive stars. The young stars thus formed will 
be roughly concentrated between the remnant cloud and the 
massive stars. In Fig. [7j the spatial distribution of YSOs in 
NGC 6823 shows a concentration of Class II sources between 
the pillar VulP12 in the north-east and the OB stars at the 
center. The pillar seems to be pointed towards the massive 
stars at the center. This suggests that the RDI mechanism 
may have triggered the formation of the YSOs seen in the 
eastern part of the cluster. In contrast, the western part does 
not contain any pillars or any such extended structures, and 
is also found to be nearly devoid of disk sources. 
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Such triggered star formation is self-sustaining in time 
and self-propagating in space, as compared to spontaneous 
cloud collapse. The excess of YSOs mainly to the east of the 
Trapezium suggests the presence of some external stimulus, 
which could be the RDI mechanism. 



5 SUMMARY 

We present results from an optical through MIR photomet- 
ric survey of the young open cluster NGC 6823 that lies in 
the Vulpecula OBI association at ~2kpc. There is signifi- 
cant differential reddening in the cluster, with Av ranging 
between ~3-20 mag. We find a ~20% fraction of Class I/II 
sources in the cluster. There is also a large population of 
photospheric Class III sources, a majority of which is likely 
the extincted background/foreground field star population. 
A large fraction of the disk sources are found to be among 
the younger group of low-mass stars in NGC 6823, at ages 
of ~l-5 Myr and masses of ~0.1-0.4Mq. The age and mass 
estimates are tentative and require cluster membership con- 
firmation. A majority of the disk sources are located in the 
eastern part of the cluster, while the western region mostly 
consists of the Class Ill/field stars. Our survey has been 
deep enough to detect the low-mass population in this clus- 
ter. NGC 6823 provides a modest density comparison to the 
Orion Nebula Cluster, which also has an extreme environ- 
ment with both a high- and a low-mass population extending 
from massive 0-type stars down to sub-stellar objects (e.g., 
Slesnick et al. 2004). The presence of young disk sources 
in NGC 6823 indicates similar star formation properties in 
the outer regions of the Galaxy as observed for young clus- 
ters in the solar neighborhood. Follow-up studies including 
multi-object spectroscopy and X-ray imaging are needed to 
deepen our understanding of this cluster. 
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Table 1. Optical, NIR, and MIR Photometry 
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Table 1. Continued 



RA (J2000)Deo (J2000) 



[3.6] 



[4.5] 



[5.8 



[S] 



19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43; 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43 
19:43; 
19:43; 



13.02 +23; 
13.10 +23; 

13.14 +23 

13.15 +23; 
13.34 +23; 

13.34 +23; 
13.47 +23; 

13.55 +23; 

13.56 +23; 
13.62 +23; 

13.73 +23; 
13.78 +23; 
13.98 +23; 
14.28 +23 
14.30 +23; 
14.45 +23 

14.51 +23 

14.52 +23 
14.62 +23 
14.80 +23; 
14.95 +23; 

14.95 +23; 

14.96 +23; 
14.96 +23; 

14.96 +23; 
15.09 +23; 

15.14 +23 
15.68 +23; 

15.74 +23; 

15.78 +23; 
15.82 +23; 
16.21 +23 
16.49 +23; 

16.49 +23; 
16.64 +23; 
16.80 +23; 
16.80 +23; 

16.87 +23; 
16.92 +23; 
16.94 +23; 

16.98 +23; 

16.99 +23; 
16.99 +23; 

17.08 +23; 
17.37 +23; 
17.37 +23; 
17.47 +23; 
17.73 +23; 

18.01 +23; 

18.09 +23 

18.50 +23 

18.51 +23 
18.51 +23 

18.53 +23 
18.89 +23 

19.02 +23 

19.03 +23 

19.15 +23 
19.30 +23 
19.49 +23 
19.55 +23 
19.60 +23 
19.60 +23 
19.72 +23 

19.92 +23 

19.93 +23 
20.03 +23 
20.20 +23 
20.24 +23 
20.49 +23 
20.70 +23 

20.88 +23 

20.97 +23 

20.98 +23 
21.06 +23 
21.13 +23 
21.15 +23 
21.19 +23 
21.28 +23 
21.28 +23 

21.47 +23 
21.68 +23 

21.79 +23 
21.82 +23 

22.48 +23 

22.57 +23 
22.57 +23 

22.62 +23 

22.63 +23 
22.77 +23 
22.85 +23 

22.97 +23 
23.00 +23 

23.03 +23 
23.17 +23 
23.24 +23 
23.33 +23 

23.35 +23 
23.39 +23 
23.57 +23 
23.57 +23 

23.98 +23 

23.99 +23 

24.04 +23 
24.04 +23 
24.06 +23 
24.12 +23 



16:53.18 

19:59.88 

15:38.09 

13:33.67 

14:34.48 

14:33.43 

19:15.78 

16:53.00 

17:58.74 

17:40.16 

21:5.40 

21:14.51 

15:58.43 

21:10.84 

17:17.41 

17:59.42 

20:15.79 

21:1.12 

16:0.77 

15:30.28 

13:56.60 

21:39.64 

16:59.88 

16:58.94 

21:41.40 

17:20.04 

17:20.83 

17:6.07 

17:9.13 

17:6.07 

21:39.89 

15:53.57 

21:29.88 

21:29.30 

17:43.40 

20:47.11 

19:17.62 

19:10.81 

18:4.28 

17:58.88 

17:51.65 

19:48.50 

19:37.16 

14:33.00 

20:25.48 

20:26.12 

18:12.35 

17:21.77 

17:44.84 

15:40.36 

20:45.49 

20:33.90 

20:44.88 

17:4.24 

18:44.46 

21:56.84 

19:24.46 

20:29.11 

18:47.56 

21:66.95 

20:33.47 

19:10.63 

19:9.88 

17:46.25 

18:0.00 

19:0.26 

17:25.40 

17:6.47 

19:19.78 

18:31.82 

19:9.26 

17:31.13 

19:1.63 

19:2.39 

21:28.73 

19:30.47 

19:57.72 

17:33.79 

14:40.13 

16:12.50 

16:25.28 

19:54.23 

17:32.60 

21:4.07 

18:24.55 

16:31.91 

16:30.94 

21:7.56 

21:7.38 

18:26.60 

18:49.82 

17:59.71 

18:9.97 

20:54.92 

19:34.21 

20:35.77 

17:43.87 

21:40.10 

21:38.56 

18:37.19 

20:8.81 

20:7.66 

21:57.24 

13:56.31 

16:26.39 

15:2.09 

19:22.80 



20.38(0.03) 
19.68(0.03) 
18.17(0.03) 

15.94(0.02) 
15.47(0.02) 



17.68(0.01) 
10.99(0.02) 

14.11(0.01) 

11.28(0.03) 
20.99(0.03) 
19.75(0.06) 
14.52(0.02) 
21.71(0.03) 
21.01(0.03) 

15.64(0.02) 

19.09(0.02) 

19.79(0.02) 

16.24(0.03) 

12.76(0.02) 

12.18(0.01) 

20.47(0.02) 
11.86(0.01) 
13.33(0.01) 
21.22(0.04) 

15.05(0.01) 
18.63(0.04) 



15.44(0.01) 
14.27(0.02) 
13.12(0.01) 



15.44(0.02) 
22.16(0.08) 
19.18(0.01) 



21.34(0.03) 



12.77(0.01) 
17.92(0.01) 



17.48(0.01) 
18.18(0.01) 
19.11(0.01) 



21.04(0.02) 



22.66(0.04) 

13.00(0.02) 
21.36(0.02) 
12.42(0.02) 



15.22(0.02) 



15.04(0.02) 



13.57(0.02) 
18.14(0.01) 



20.22(0.02) 
21.88(0.04) 



16.30(0.02) 
18.47(0.02) 



16.85(0.01) 



20.47(0.02) 
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24 Riaz et al. 



Table 2. MIPS Photometry 


RA (J2000) 


Dec (J2000) 


[24]^] 


SED Class 


1 Q-ZL9- '^7 






11 






5.07 


jj 


1 Q-A'\-C\A AO 

± y . . . 






jj 


19:43:06.46 


+23:14:50.96 


7.31 


II 


19:43:08.33 


+23:19:13.01 


3.49 


II 


19:43:09.97 


+23:18:01.08 


7.3 


II 


19:43:11.58 


+23:17:32.42 


7.37 


II 


19:43:11.58 


+23:18:20.45 


7.37 


II 


19:43:12.47 


+23:18:17.03 


8.03 


III 


19:43:13.47 


+23:19:15.78 


6.97 


II 


19:43:13.56 


+23:17:58.74 


7.76 


II 


19:43:14.95 


+23:13:56.6 


6.99 


II 


19:43:20.98 


+23:19:02.39 


3.73 


II 


19:43:24.42 


+23:20:12.16 


7.23 


II 


19:43:24.80 


+23:18:53.35 


7.57 


II 


19:43:30.06 


+23:15:05.65 


6.75 


II 



The uncertainty on the 24|xm photometry is estimated to be <0.2 mag (Billot et al. 
2010). 
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Table 3. IPHAS Photometry 



RA (J2000) Dec (J2000) 



Ha 



SED Clasf " 



19:42 


47 


28 


+23:19:35.00 


17 


78 


(0 


008) 


15 


47 


(0 


005) 


17 


11 


(0 


014) 


III 


19:42 


47 


39 


+23:18:24.98 


20 


15 


(0 


045) 


17 


84 


(0 


023) 


19 


66 


(0 


119) 


III 


19:42 


47 


49 


+23:18:54.14 


19 


05 


(0 


019) 


17 


47 


(0 


018) 


18 


53 


(0 


045) 


III 


19:42 
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61 


+23:17:1.32 


17 


80 


(0 


008) 


15 


67 


(0 


006) 


17 


14 


(0 


013) 


III 


19:42 
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64 


+23:19:51.82 


20 


26 


(0 


049) 


17 


97 


(0 


025) 


19 


41 


(0 


096) 


III 


19:42 


47 


65 


+23:17:9.24 


21 


01 


(0 


089) 


18 


52 


(0 


037) 


20 


21 


(0 


157) 


III 


19:42 


48 


45 


+23:17:30.01 


19 


94 


(0 


036) 


17 


69 


(0 


020) 


19 


28 


(0 


069) 


III 


19:42 


48 


49 


+23:17:30.01 


19 


87 


(0 


036) 
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65 
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020) 


19 


16 


(0 


077) 


III 
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60 
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73 
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014) 
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02 
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012) 
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19 


(0 


028) 


III 
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90 
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20 


69 


(0 


067) 


18 


31 


(0 


032) 


19 


86 
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115) 
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19:42 


49 


39 


+23:17:28.28 


18 


62 
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013) 


16 


46 


(0 


009) 


17 


99 
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024) 
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19:42 


49 


45 


+23:14:33.68 


14 


94 
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002) 


14 


31 
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003) 


14 


65 


(0 


003) 
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49 


63 


+23:17:0.49 


20 


11 


(0 


041) 


17 


15 


(0 


014) 


19 


30 


(0 


071) 


II 


19:42 


49 


64 


+23:18:40.86 


18 


91 


(0 


016) 


17 


85 


(0 


022) 


18 


34 


(0 


032) 


III 


19:42 


49 


85 


+23:15:58.03 


20 


58 


(0 


065) 


18 


10 


(0 


028) 


20 


03 


(0 


166) 


III 


19:42 
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17 


+23:15:38.99 


13 


48 
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001) 
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06 
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002) 


13 


25 
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002) 
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17 
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13 


65 
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001) 
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10 
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39 
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002) 
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20 
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04 
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41 
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003) 
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71 
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003) 


III 


19:42 
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08 
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13 


37 
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001) 
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95 
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001) 


13 


19 
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002) 
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15 
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92 
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016) 
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54 
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018) 
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14 
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027) 


II 


19:42 
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17 
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21 


16 


(0 


102) 


17 


33 


(0 


015) 


19 


63 


(0 


094) 


III * 


19:42 


52 


30 


+23:19:56.53 


17 


43 


(0 


006) 


16 


33 


(0 


008) 


16 


85 


(0 


Oil) 


III 


19:42 


52 


68 


+23:18:20.59 


15 


88 


(0 


003) 


14 


98 


(0 


004) 


15 


21 


(0 


004) 


II * 


19:42 


52 


86 


+23:15:17.06 


20 


28 


(0 


048) 


18 


76 


(0 


045) 


19 


26 


(0 


068) 


II * 


19:42 


52 


98 


+23:14:6.40 


14 


74 


(0 


002) 


14 


15 


(0 


003) 


14 


44 


(0 


003) 


III 


19:42 


53 


18 


+23:13:51.60 


21 


30 
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115) 


19 


03 
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057) 
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59 
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091) 
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53 


25 
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69 
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002) 
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96 
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29 
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83 
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59 


(0 


002) 
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(0 


004) 
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(0 


004) 
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53 


90 


+23:13:51.78 


21 


47 


(0 


135) 


18 


93 


(0 
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19 


92 


(0 
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45 


(0 


Oil) 


16 


07 


(0 


007) 


17 
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63 
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45 
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56 
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109) 
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97 
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90 
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002) 
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003) 
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53 
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003) 


III 


19:42 


55 


52 


+23:13:18.05 


20 


52 


(0 


058) 


17 


87 


(0 


023) 


19 


76 


(0 


105) 


III 


19:42 


55 


75 


+23:18:25.52 
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19:43:13.00 


1 OO "IT CrO T/1 

+23:17:58.74 


15 


60 


(0 


002) 


14 


64 


(0 


003) 


15 


20 


(0 


004) 


II 


"in /IO IO r'o 


1 OO "IT /in -\ r- 

+23:17:40.16 


15 


15 


(0 


002) 


14 


27 


(0 


003) 


14 


78 


(0 


003) 


III 


in /IO 1/1 no 

19:43:14.28 


1 OO 01 "ino/1 

+23:21:10.84 


17 


44 


(0 


006) 


16 


31 


(0 


008) 


16 


96 


(0 


013) 


III 


in./io.i/i CI 

19:4o;14.ol 


1 OO on 1 c Tn 

+23:20:15.79 


13 


85 


(0 


001) 


13 


05 


(0 


002) 


13 


52 


(0 


002) 


III 


1 n. Q . 1 /( on 

19:4o:14.o(J 


1 00.1 CT.on oo 

+23:15:30.28 


19 


76 


(0 


032) 


18 


17 


Jo 


030) 


18 


55 


(0 


045) 


II * 


in. /(Q.I nrr 

19:4cS:14.9o 


1 oo.io.rrc ^jn 

+23:13:56.60 


18 


19 


(0 


010) 


16 


67 


Jo 


010) 


17 


58 


(0 


021) 


II 


in. /(O.I /( ncr 

19:4o:14.9o 


1 oo.m .on a A 

H-zo:zi:oy.D4 


14 


34 


(0 


001) 


13 


74 


(0 


002) 


14 


01 


(0 


002) 


III 


1 n. /( o . 1 A r\c 

19:4o:14.yD 


1 oo -\ r- trn oo 

+23:16:59.88 


20 


55 


(0 


063) 


18 


57 


Jo 


040) 


19 


78 


(0 


133) 


II 


1 n . /( o . 1 c nn 

19:43:15. u9 


1 OO "1 T on n A 

+23:17:20.04 


15 


28 


(0 


002) 


14 


38 


(0 


003) 


14 


73 


(0 


003) 


III 


"in /I O 1 CT TO 

19:43:15.78 


1 O O "1 T 1^ nT 

+23:17:6.07 


18 


81 


(0 


015) 


17 


52 


(0 


018) 


18 


29 


(0 


037) 


III 


"in /IO i(^ 0"i 

ly;4o:iD.zi 


1 OO "icr cro ctt 

+23:15:53.57 


15 


84 


(0 


003) 


14 


88 


(0 


004) 


15 


28 


(0 


005) 


III 


"in /I o 1 on 

19:43:10. 8U 


1 OO "1 n "1 T (^o 

+23:19:17.62 


18 


17 


(0 


010) 


16 


86 


(0 


012) 


17 


64 


(0 


022) 


III 


in /I o 1 no 

19:43:10.98 


1 OO IT CTI 

+23:17:51.65 


19 


63 


(0 


029) 


18 


08 


(0 


028) 


18 


87 


(0 


060) 


III 


1 n. /(O.I nn 

19:43:16.99 


1 oo.in.OT Id 

+23:19:37.16 


13 


30 


(0 


001) 


12 


72 


Jo 


001) 


13 


02 


(0 


001) 


III 


in./(o.i'7 no 

19:4cS:17.Uo 


1 00.1 /i.Oo nn 

+2c>:14:Jo.UU 


20 


29 


(0 


051) 


18 


26 


Jo 


031) 


18 


92 


(0 


064) 


II * 


iy:4o:i7.o7 




1 A 

14 


71 


(0 


002 j 


13 


80 


Jo 


002) 


1 A 

14 


15 


(0 


003) 


TTT 
ill 


in./(Q.i'7 /(*? 

19:4d:17.47 


1 00.10.10 OCT 

+2c>:lo:12.o5 


14 


39 


(0 


001) 


13 


28 


(0 


002) 


13 


90 


(0 


002) 


III 


iy:4o:lo.Ul 


1 OQ . 1 T. A A Q A 

-\-Zo.L ( :44.o4 


18 


18 


(0 


010) 


16 


87 


Jo 


012) 


17 


62 


(0 


022) 


TT 
11 


in /I o 1 o crn 

iy:4v>:lo.o(J 


1 oo on ^ r ^c\ 


15 


11 


(0 


002) 


14 


28 


(0 


003) 


14 


78 


(0 


004) 


III 


in./io."io c^Q 
iy;4o:io.oo 


1 O Q . 1 T. /I O /I 


14 


24 


(0 


001) 


13 


i 3 


Jo 


002) 


13 


96 


(0 


002) 


TTT 
lii 


in /I o 1 o on 

19:43:18.89 


1 oo "IO A A A r- 

+zo:1o:44.4d 


13 


15 


(0 


001) 


12 


65 


(0 


001) 


12 


87 


(0 


001) 


III 


19:43:19.15 


1 OQ.on.on 1 1 
H-zo:zU:zy. li 


15 


11 


(0 


002) 


14 


10 


Jo 


003) 


14 


71 


(0 


003) 


TTT 
ill 


19:43:19.49 


1 00.01 .T^£i n CT 

+23:21:56.95 


18 


43 


(0 


012) 


17 


53 


(0 


018) 


17 


93 


(0 


028) 


III 


19:43:2U.U3 


+23:17:25.40 


19 


88 


(0 


036) 


16 


76 


Jo 


Oil) 


19 


15 


(0 


078) 


TTT 
ill 


19:43:20.49 


+23:18:31.82 


17 


59 


(0 


007) 


16 


16 


JO 


008) 


17 


16 


(0 


015) 


III 


19:43:20.90 


I OQ.in.o on 
+23:19:2.39 


16 


99 


(0 


005) 


15 


71 


Jo 


006) 


16 


09 


(0 


007) 


TT * 
11 


19:43:21.06 


+23:21:28.73 


18 


81 


(0 


015) 


17 


64 


(0 


020) 


18 


03 


(0 


030) 


II * 


19:43:21.19 


+23:17:33.79 


19 


68 


(0 


030) 


17 


99 


(0 


026) 


18 


90 


(0 


061) 


II 


19:43:21.28 


+23:16:12.50 


21 


30 


(0 


123) 


18 


76 


(0 


046) 


19 


47 


(0 


103) 


III * 


19:43:21.47 


+23:16:25.28 


21 


50 


(0 


146) 


19 


10 


(0 


063) 


19 


34 


(0 


090) 


III * 


19:43:21.79 


+23:17:32.60 


13 


08 


(0 


001) 


12 


65 


(0 


001) 


12 


85 


(0 


001) 


III 


19:43:21.82 


+23:21:4.07 


19 


76 


(0 


035) 


18 


76 


(0 


053) 


19 


72 


(0 


125) 


II 


19:43:22.62 


+23:21:7.56 


14 


89 


(0 


002) 


14 


15 


(0 


003) 


14 


48 


(0 


003) 


III 


19:43:22.85 


+23:18:49.82 


14 


94 


(0 


002) 


14 


36 


(0 


003) 


14 


61 


(0 


003) 


III 


19:43:23.00 


+23:18:9.97 


13 


70 


(0 


001) 


13 


38 


(0 


002) 


13 


56 


(0 


002) 


III 


19:43:23.03 


+23:20:54.92 


17 


84 


(0 


009) 


16 


57 


(0 


Oil) 


17 


21 


(0 


016) 


III 


19:43:23.33 


+23:17:43.87 


20 


58 


(0 


065) 


17 


53 


(0 


018) 


19 


65 


(0 


121) 


III 


19:43:23.39 


+23:21:38.56 


16 


21 


(0 


003) 


15 


31 


(0 


005) 


15 


73 


(0 


006) 


II 
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Table 3. Continued 



RA (J 2000) 


Dec (J 2000) 




r 


r ] 




r 








[Ha] 


abjD Class 


19:43:23.57 


+23:18:37.19 


18 


27 


(0.011) 


17 


09 


(0 


013) 


17 


49 


(0.019) 


II * 


19:43:23.98 


+23:20:7.66 


16 


58 


(0.004) 


15 


50 


(0 


005) 


16 


05 


(0.007) 


III 


19:43:24.12 


+23:19:12.11 


19 


09 


(0.019) 


18 


27 


(0 


032) 


18 


62 


(0.048) 


III 


19:43:24.13 


+23:19:23.59 


20 


20 


(0.047) 


18 


50 


(0 


038) 


19 


46 


(0.100) 


II 


19:43:24.16 


+23:17:24.97 


19 


30 


(0.022) 


17 


93 


(0 


024) 


18 


40 


(0.041) 


II * 


19:43:24.42 


+23:20:12.16 


17 


27 


(0.006) 


16 


30 


(0 


008) 


16 


64 


(0.011) 


II 


19:43:24.52 


+23:21:13.03 


13 


05 


(0.001) 


12 


67 


(0 


001) 


12 


76 


(0.001) 


III 


19:43:24.71 


+23:19:11.68 


15 


17 


(0.002) 


14 


57 


(0 


003) 


14 


84 


(0.004) 


III 


19:43:24.80 


+23:18:53.35 


18 


98 


(0.018) 


17 


58 


(0 


019) 


17 


39 


(0.018) 


II * 


19:43:25.10 


+23:18:6.41 


20 


08 


(0.042) 


17 


48 


(0 


018) 


19 


43 


(0.097) 


III 


19:43:25.18 


+23:15:56.34 


17 


16 


(0.005) 


15 


97 


(0 


007) 


16 


58 


(0.010) 


III 


19:43:25.23 


+23:19:1.52 


19 


52 


(0.027) 


18 


18 


(0 


030) 


18 


74 


(0.053) 


III 


19:43:25.45 


+23:20:17.34 


16 


51 


(0.004) 


14 


77 


(0 


004) 


15 


84 


(0.007) 


III 


19:43:25.61 


+23:19:22.30 


15 


26 


(0.002) 


14 


38 


(0 


003) 


14 


80 


(0.003) 


III 


19:43:25.93 


+23:22:7.07 


18 


79 


(0.015) 


17 


15 


(0 


014) 


18 


30 


(0.037) 


III 


19:43:26.19 


+23:20:31.45 


13 


67 


(0.001) 


12 


78 


(0 


001) 


13 


20 


(0.002) 


III 


19:43:26.30 


+23:21:18.90 


20 


52 


(0.067) 


19 


24 


(0 


079) 


19 


28 


(0.085) 


I 


19:43:26.44 


+23:15:28.66 


18 


98 


(0.018) 


17 


72 


(0 


021) 


17 


94 


(0.027) 


II * 


19:43:27.01 


+23:19:44.44 


18 


40 


(0.012) 


17 


10 


(0 


014) 


17 


38 


(0.018) 


II * 


19:43:27.32 


+23:18:41.00 


18 


51 


(0.012) 


16 


03 


(0 


007) 


17 


66 


(0.022) 


III 


19:43:28.20 


+23:21:34.38 


17 


16 


(0.006) 


15 


57 


(0 


006) 


16 


57 


(0.010) 


III 


19:43:28.26 


+23:18:48.20 


20 


66 


(0.069) 


18 


87 


(0 


052) 


19 


14 


(0.075) 


II * 


19:43:28.76 


+23:18:25.56 


16 


36 


(0.003) 


15 


31 


(0 


005) 


15 


y ( 


,U.UU ( ) 


III 


19:43:29.05 


+23:21:23.22 


20 


33 


(0.052) 


18 


33 


(0 


033) 


18 


01 


(0.029) 


III * 


19:43:30.03 


+23:20:53.02 


14 


36 


(0.001) 


13 


57 


(0 


002) 


13 


97 


(0.002) 


III 


19:43:30.06 


+23:15:5.65 


18 


47 


(0.012) 


17 


36 


(0 


016) 


17 


33 


(0.017) 


II * 


19:43:30.19 


+23:21:18.22 


14 


20 


(0.001) 


13 


47 


(0 


002) 


13 


80 


(0.002) 


III 


19:43:30.58 


+23:22:31.94 


19 


89 


(0.039) 


17 


31 


(0 


017) 


18 


70 


(0.052) 


III 


19:43:30.70 


+23:15:54.22 


21 


38 


(0.131) 


19 


34 


(0 


075) 


19 


45 


(0.101) 


III * 


19:43:30.74 


+23:17:50.06 


16 


14 


(0.003) 


15 


38 


(0 


005) 


15 


73 


(0.006) 


III 


19:43:31.18 


+23:19:23.99 


19 


96 


(0.038) 


17 


57 


(0 


019) 


19 


62 


(0.115) 


III 


19:43:31.20 


+23:19:28.56 


15 


22 


(0.002) 


14 


67 


(0 


003) 


14 


91 


(0.004) 


III 


19:43:31.87 


+23:16:50.95 


15 


97 


(0.003) 


15 


37 


(0 


005) 


15 


61 


(0.005) 


II 


19:43:31.96 


+23:16:55.60 


20 


25 


(0.049) 


18 


18 


(0 


029) 


19 


68 


(0.124) 


III 


19:43:32.20 


+23:20:21.41 


17 


75 


(0.008) 


16 


73 


(0 


Oil) 


17 


04 


(0.014) 


II * 


19:43:32.21 


+23:17:36.89 


19 


34 


(0.023) 


17 


37 


(0 


016) 


18 


87 


(0.061) 


III 



